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PACS 74 . 70 . Xa - Pnictides and chalcogenides 

PACS 74 . 25 . Bt - Thermodynamic properties 

PACS 65 . 40 . Ba - Heat capacity 

PACS 74 . 20 . Rp - Pairing symmetries (other than s-wave) 

Abstract. - An extensive calorimetric study of the normal- and superconducting-state proper- 
ties of Ba(Fei_ 2 ,Co 2 ,)2As2 is presented for < x < 0.2. The normal-state Sommerfeld coefficient 
increases (decreases) with Co doping for x < 0.06 (x > 0.06), which illustrates the strong competi- 
tion between magnetism and superconductivity to monopolize the Fermi surface in the underdoped 
region and the filling of the hole bands for overdoped Ba(Fei_ I Co a; )2As2. All superconducting 
samples exhibit a residual electronic density of states of unknown origin in the zero-temperature 
limit, which is minimal at optimal doping but increases to the normal-state value in the strongly 
under- and over-doped regions. The remaining specific heat in the superconducting state is well 
described using a two-band model with isotropic s-wave superconducting gaps. 



Introduction. — Despite the fact that the theoretical 
background has been available since the late 50's with the 
pioneering papers of Suhl et al. [1] and Moskalenko [2], 
multiband superconductivity (MBSC) emerged as an 
unanimously accepted phenomenon only after the discov- 
ery of the MgB2 superconductor in 2001 [3]. Rapidly, 
calorimetric signatures, like the excess specific heat ob- 
served at low temperature (with respect to the single-band 
BCS curve) [I], the initial rapid rise of the mixed-state 
heat capacity with magnetic field 5 , and the anomalous 
positive curvature of the upper critical field 6 provided 
the most convincing evidence of its existence and are now 
textbook hallmarks of the existence of two gaps. Later, 
significant interband contributions to the Eliashberg func- 
tion, reminiscent of MBSC, were experimentally detected 
using tunneling experiments |49j . Since then, the occur- 
rence of MBSC has been discussed for many different com- 
pounds including heavy fermions [7] , cobaltates [8] , chalco- 
genides [9], A15 compounds [50], and the recently discov- 
ered iron-pnictide family. The aforementioned character- 
istic signatures of MBSC arc less pronounced in iron pnic- 
tide superconductors since the interband coupling, pro- 



moted by the strong Fermi-surface nesting, is more im- 
portant than in MgB2 leading to a smaller gap anisotropy. 
Although, the existence of several energy gaps appears to 
be accepted, the precise symmetry of the superconducting 
order parameter and its evolution with doping remain con- 
troversial. Specific-heat [10l[TT] and ARPES [H1[T3] stud- 
ies agree on isotropic s-wave gaps in Ba(Fei_ :r Co a ;)2As2 
and Bai_ £C K; C Fe2As2, but differ on the gap amplitudes. On 
the other hand, Raman scattering [14] and heat transport 
measurements [15] . which are both interpreted in terms 
of accidental nodes (and/or gap minima), show contradic- 
tory behavior with doping. 

In this Letter, we present an extensive thermodynamic 
study of the evolution of the normal- and superconduct- 
ing state properties of Ba(Fei_ a; Co a ;)2As2 single crystals 
with Co concentration in the range < x < 0.2. We 
find that the normal-state Sommerfeld coefficient (7„) 
first increases with Co doping up to x = 0.06 (where T c 
is maximal), which is related to the gradual suppression 
of the spin-density-wave (SDW) state; j„ then decreases 
for x > 0.06, which we attribute to a gradual filling of 
the hole bands with increasing Co content. This illus- 
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trates the strong competition between magnetism and su- 
perconductivity to monopolize the Fermi surface (FS) in 
Ba(Fei_ K Coa;)2As2. Further, we find that all supercon- 
ducting samples exhibit a residual electronic density of 
states (7 r ) in the zero-temperature limit, which is mini- 
mal at optimal doping but increases up to the normal-state 
value in the strongly under- and overdoped states. We 
show that this residual term is intrinsic and not due to, 
e.g., residual flux or oxidation products in the crystals. 
The observed values of j r imply that the superconduct- 
ing fraction strongly depends on doping and is only close 
to unity near optimal doping. We compare this behavior 
with that of the cuprates, which show a similar tendency. 
The superconducting-state specific heat is well described 
using a two-band model with isotropic s-wave gaps. The 
smaller gap ratio Ai(0)/kaT c is nearly independent of 
doping, whereas the large gap ratio A 2 (0)/ksT c shows a 
tendency to strong coupling at optimal doping. 

Experimental details. — Single crystals, with typi- 
cal mass of 5-10 mg, were grown from self- flux using pre- 
reacted FeAs and CoAs powders mixed with Ba. Details 
of the crystal synthesis are described elsewhere [TUIfTo] . 
The precise composition of all samples was determined by 
energy dispersive x-ray spectroscopy and complemented 
by 4-circle diffractometry. As reported by other groups, 
our characterization shows that the effective Co amount 
in the sample represents only 75 % of the nominal flux 
Co content. The specific heat was measured by relax- 
ation calorimetry down to 0.4 K, with the 3 He-PPMS from 
Quantum Design. We have observed a slow degradation 
of the samples with time (over several months), which 
resulted in the progressive growth of a low-temperature 
Schottky anomaly. Therefore, all the measurements pre- 
sented in this article were performed on crystals which 
were cooled down and measured within less than 10 hours 
after opening the Si02 ampoule. 

Specific-heat measurements. — Figures [TJa)-(b) 
display the specific heat plotted as C/T as a function 
of T 2 (above 0.4 K) of non-superconducting underdoped 
(x < 0.04) and overdoped (x > 0.14) single crystals be- 
low 6 K. In contrast to other reports, we do not observe 
any Schottky-like contribution, revealing the high quality 
of these samples. The heat capacity accurately follows 
the standard behavior, C/T=7 ra +B3T 2 , where j n and B 3 
are the normal-state Sommcrfeld and Debye terms, re- 
spectively. The Debye temperature changes by less than 
3 % between the low-doping orthorhombic and the over- 
doped tetragonal phases confirming that the substitution 
of Fe by Co barely affects the phonon spectrum, in ac- 
cord with neutron scattering measurements [17] . Figures 
[TJc)-(j) show the electronic specific-heat contribution (C e ) 
of superconducting crystals with 0.04 < x < 0.12, ob- 
tained by subtracting the lattice specific heat from the 
raw data following the procedure given in Ref. [lOj , For 
overdoped (underdoped) single crystals, we used the lat- 
tice contribution derived from the x = 0.153 (x = 0.03) 
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Fig. 1: Specific heat of underdoped (yellow symbols), opti- 
mally doped (magenta symbols) and overdoped (blue sym- 
bols) Ba(Fei- i rCo : r)2As2 single crystals, (a)-(b) C/T of non- 
superconducting samples as a function of T 2 . Lines are fit 
to the Debye model with j n , B3, and Qd given in mj mol -1 
K~ , mj mol - K~ , and K, respectively, (c)-(j) Electron spe- 
cific heat C e /T of several superconducting samples. Dotted 
and dashed lines represent 7n and j r , respectively (numerical 
values are given). Solid lines are fits to the two-band a-model. 



sample. All data comply with entropy conservation, i.e., 
J Tc {Ce/T) dT = 7n T c . Except for the x = 0.04 and x 
= 0.12 curves, the superconducting transitions at T c are 
remarkably sharp, indicating little inhomogeneity in the 
crystals. Nevertheless, our results show that C e /T does 
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not extrapolate to zero at T = 0, but exhibits a residual 
normal-state-like contribution, r y r . In the following, we 
discuss in detail the doping evolution of (i) the normal- 
state electronic density of states (EDOS) and the resulting 
phase diagram, (ii) the superconducting gaps derived from 
fitting the specific heat with a two-band model (solid lines 
in Fig. QJc)-(j)) and (hi) the residual electronic density of 
states (RDOS). 

(x,T) phase diagram and normal-state EDOS. 

Fig. EJa) shows the phase diagram derived from our 
specific-heat data together with thermal-expansion mea- 
surements performed on the same samples (Ref. [18]). 
Our phase diagram is very similar to previous reports, 
except that superconductivity disappears at a lower Co 
concentration (x « 0.125) than previously reported from 
transport and magnetization measurements [T9"M22| . We 
also observe that superconductivity emerges already in the 
magnetic phase, at a concentration x ss 0.04, and T c in- 
creases to 24.2 K at x w 0.575. This increase is accompa- 
nied by a rise of 7„ from 5.6 mJ mol -1 K~ 2 at x = to 
about 22 mJ mol -1 K~ 2 close to the optimal concentra- 
tion (see Fig. HJb)), which we attribute to the closing of 
the SDW gap with Co doping. Our data show that the de- 
crease of 7 n with increasing x is linear for x > 0.0575 and 
extrapolates for x = to a value of « 27 mJ mol -1 K~ 2 
(dashed line in Fig. EJb)), which we interpret as the high- 
temperature estimate of the Sommerfeld coefficient for T 
> T^r, i.e., in the paramagnetic phase of BaFe2As2- This 
indicates that the EDOS is reduced by a factor of six due 
to the FS reconstruction caused by the SDW. An increase 
of 7 n implies that a larger fraction of the Fermi surface 
becomes available for superconductivity, which naturally 
explains the increase of T c and AC/7„T C with Co dop- 
ing in the underdoped region. In fact, it was recently 
shown that SDW and superconducting states coexist and 
compete for the same FS in this region of the phase di- 
agram [23l[24]. On the other hand, for x > 0.575, both 
7„ and T c decrease as the system moves away from the 
magnetic instability. In Al- and C- doped MgB 2 [25], the 
reduction of j n with electron doping is well understood 
by the filling of the 2D a hole bands and this interpre- 
tation may also hold for overdoped Ba(Fei_- c Co; C )2As2. 
Indeed, ARPES measurements [26] have shown that x = 
0.10 superconducting samples have hole bands only par- 
tially filled, whereas the same bands are completely shifted 
below the Fermi level E^ when the Co content is increased 
to x = 0.15 where T c is zero. Therefore, the collapse of 
superconductivity around x = 0.13 is likely related to the 
disappearance of the interband nesting condition that con- 
nect hole and electron bands. Moreover, in doped MgB 2 , 
the reduction of the EDOS at Ep also leads to a weakening 
of the electron-phonon coupling resulting in the disappear- 
ance of superconductivity [25] . 

Two-gap analysis. — Although, in principle, three 
or more bands should be used to describe the thermody- 
namics of the Fe pnictidcs, simple two-band models have 



been successfully applied to describe a variety of exper- 
imental data. In the following, we therefore use a two- 
gap model to analyze the heat-capacity data, which has 
the advantage of a minimal number of fit parameters. As 
shown in Figs. [2c)-(j), C e is very well represented, for 
all superconducting samples, by two-gap fits similar to 
those made for MgB2 27 . In these fits, C e is the sum 
of three contributions: the residual term, 7 r T, and the 
contributions of the two gapped components of the super- 
conducting condensate derived from the two bands, i = 1, 
2. The temperature dependences of the superconducting- 
state contributions are those given by the a- model [2"9"j . 
in which the weak-coupling BCS temperature dependence 
of the gaps is assumed, but the T = value, Aj(0), is an 
adjustable parameter represented by on = Ai(0)/keT c (in 
the BCS weak-coupling limit a(0) = 1.764). The ampli- 
tudes of these contributions are proportional to ji, with 
the following constraint: 71 + 72 = 7n - lr- Thus, the fits 
introduce three adjustable parameters, a%, ot<2, and 71/72. 

Fig. Wis) shows the specific-heat ratio, AC/7„T C , as 
a function of Co content. At the optimal concentra- 
tion, it reaches a value of about 1.6, which slightly ex- 
ceeds the weak-coupling BCS value, but is significantly 
smaller than 2.5 reported for Bao.6sKo.32Fe2As2 [IT]. The 
evolution of AC/7„T C and 7 r with doping is strongly 
correlated, i.e. the largest (smallest) AC/7„T C occurs 
for the smallest (largest) j r . This indicates that the 
strong doping dependence of AC/7„T C is due primar- 
ily to the presence of j r . Indeed, after normalizing 
AC to the superconducting fraction (7 n -7 r ) instead of 
7„, the ratio AC/(7„ — 7 r )T c exhibits a weaker dop- 
ing dependence than AC/7„T C , which likely reflects the 
evolution of A2(0)/kgT c with x. We note that recent 
specific-heat measurements by Bud'ko et al. [30] have 
reported that AC/T C ex T 2 for Ba(Fei_ 2 ,Co a; )2As2 and 
Ba(Fei_ a; Ni a ;)2As2, which was subsequently interpreted by 
Zaanen j31j in terms of a quantum critical metal under- 
going a pairing instability. We find that the same scaling 
relation holds, but with an exponent closer to ~ 1.5 (not 
shown), and it is not clear whether this theoretical ap- 
proach can explain the present data. In particular, any 
valid theory needs to explain that AC is reduced because 
of the large 7,. , and also needs to consider that in a two- 
band superconductor, the size of the jump is governed by 
both the gap anisotropy, A2/A1, and the ratio of the dif- 
ferent band densities of states, 72/71- 

Figures EJd)-(e) summarize the fit parameters derived 
from the isotropic two-band analysis. Except for the ex- 
tremely underdoped and overdoped samples, the gap val- 
ues (see Fig. [2Jd)) conform to the theoretical constraints 
that, for a weakly coupled two-band superconductor [32] . 
one gap must be larger than the BCS value and the other 
smaller. The largest gap values are obtained at optimal 
doping with a small gap 2Ai(0)/keT c « 1.9 and a major 
one with 2A2(0)/ksT c « 5. Recently, the same analy- 
sis was successfully applied to the heat capacity of the 
close to optimally doped Bao.6sKo.32Fe2As2 by Popovich 
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et al. [TT]. They found a comparable 2Ai(0)/k B T c ra 2.2, 
whereas the larger ratio, 2A 2 (0)/ksT c rj 6.6, is substan- 
tially higher, which was interpreted as a possible sign of 
strong coupling. 

FigMe) illustrates that 72 is always greater than 7! in- 
dicating that the major gap develops around the Fermi 
surface sheet that exhibits the largest EDOS, while 
in a two-band model it is theoretically expected that 
72 /71 oc -\/Ai/A2 in the limit of pure interband coupling 
within the weak-coupling limit [33j . These findings im- 
ply that either intraband couplings are not negligible or 
that more than two bands are involved in Co-doped 122 
pnictides. As pointed out by Benfatto et al. [34], a min- 
imal model of three bands is indeed required to describe 
the superconducting-state properties, in agreement with 
ARPES measurements in Bao.6Ko.4Fe2As2 [13] which have 
revealed the existence of three energy gaps. Our data 
can also be fitted with three gaps, but this would necessi- 
tate five fitting parameters. Our present analysis can be 
qualitatively reconciled with photoemission data by as- 
suming that a small gap (2Ai(0)/k_BT c rj 1.9) exists on 
the quasi-3D hole band (a band), while two large gaps 
of equal amplitude (2A 2 (0)/kBT c w 5) develop around 
the two strongly nested sheets, i.e., the hole-like sheets 
((3 bands) and the electron- like bands (7 bands). Recent 
first-principles calculations by Nakai et al. [35j support 
this assumption. In this context, 71 represents the density 
of states of the a band while 72 consists of contributions 
from the 2D j3 and 7 sheets. In the following, we discuss 
our data using this scheme. 

In Fig(3] we plot both gap amplitudes as a function of 
T c together with values for Bao.6sKo.32Fe2As2 [TT]. We 
find that (i) the small gap increases linearly with T c , (ii) 
the large gap increases stronger than linearly with T c , 
and (iii) the same behavior occurs for both under- and 
overdoped samples. Additionally, the K-doped data fit 
nicely onto the same curves. In the limit of T c — > 0, the 
larger gap A 2 (0) converges to the BCS value, Abcs(0) 
= 1.764ksT c , which has previously been reported also 
for Mgi_ a; Al a ;B2 |36J. In the latter system, the coupling 
strength and the interband scattering between the 2D a 
bands (which exhibit the larger gap) and the 3D ir bands 
are weak because of their difference in character. By anal- 
ogy, the same argument can also explain the difference 
in gap amplitude of the quasi-3D a band and the almost 
2D j3 and 7 bands in Ba(Fei_ a; Co : i;)2As2. Figure [3] shows 
a tendency for strong-coupling effects for the large gap, 
i.e., the compounds with the highest T c , especially for the 
K-doped superconductor. Recent calculations [37^ have 
shown that the Lindhard function, which provides a mea- 
sure of the degree of nesting, is much higher for K- than 
for Co-doped Bal22, and this strong coupling may thus 
be related to the degree of nesting. Taken together, our 
data and analysis are fully consistent with superconduc- 
tivity driven by interband coupling between the nested /3 
and 7 bands. 



Residual density of states, 7,.. — We now focus our 
discussion on the RDOS shown in Fig. [DJb). Sizeable j r 
values appear to be a general feature of specific-heat mea- 
surements on Co-doped 122 iron pnictides [TU1,|2"51I35] . and 
the first question to be addressed is: Is this an intrinsic or 
an extrinsic (e.g., residual flux, a second phase, or gener- 
ally poor sample quality) effect? The systematic doping 
dependence of 7 r shown in Fig. [2fb) strongly suggests an 
intrinsic origin. This is supported by recent, albeit sparse 
data of Mu et al. [28] (see Fig. Efb)), which are practi- 
cally identical in value to ours. Intrinsic RDOS is usually 
not observed in classical BCS superconductors, not even 
for dirty MgB2 [39-41]. On the other hand, many uncon- 
ventional superconductors appear to exhibit such a term. 
For example, even the best YBCO samples near optimal 
doping have -fr/ln of rougly 10 % [32], the origin of which 
remains unclear. Remarkably, an anticorrelated behavior 
between j r and 7„ similar to our results has been reported 
by Koike et al. [43] in L^-zSr^CuO^ which was argued to 
result from microscopic phase separation into SC and non- 
SC metallic regions. This could also explain our results for 
Ba(Fei_a;Coa;)2As2. A very similar explanation has been 
used in the context of penetration-depth measurements, 
where the term 'Swiss cheese' was used to describe the 
mixture of normal and superconducting regions [44j . 

On the other hand, j r could also be explained by disor- 
der, which is likely to be particularly important in uncon- 
ventional superconductors with a sign-reversing order pa- 
rameter. Recent NMR measurements [45] have shown that 
cobalt substitution does not induce any local moment, 
hence Co can be considered as a non-magnetic scatterer. 
Moreover, it was argued that the observed power-law be- 
havior of the penetration depth and the NMR spin-lattice 
relaxation time can be reconciled with an s-wave order 
paramater only if the scattering by non-magnetic impu- 
rities is taken into account [46l|47]. As shown by Onari 
and Kontani [48], interband scattering by non- magnetic 
impurities causes only a moderate drop of T c for the s++ 
state, while for a sign-reversing order parameter, like the 

s^ state, strong pair-breaking effects, including a drastic 

suppression of T c and the emergence of prominent in-gap 
states, are expected. Thus, if j r is due to disorder, then 

the order parameter is likely to be s_| However, it is not 

clear how disorder alone (due to Co doping) can explain 
the minimum of 7,. at optimal doping, since naively one 
would expect 7,. to continually increase with x. 

In summary, from an extensive thermodynamic study 
of Ba(Fei_ a; Co a ;)2As2 high-quality single crystals we find 
that the normal-state Sommerfeld coefficient first in- 
creases with Co doping up to x = 0.06 (where T c is max- 
imal) and then decreases for x > 0.06. The maximum of 
7„ at x = 0.06 arises due to the gradual suppression of the 
spin-density wave state (SDW) for x < 0.06 and to a con- 
tinuous filling of the hole bands for x > 0.06 with increas- 
ing Co content. Our results, thus, illustrate the strong 
competition between magnetism and superconductivity to 
monopolize the Fermi surface in Ba(Fei_j;Co x )2As2. All 
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our superconducting samples exhibit a residual electronic 
density of states in the zero-temperature limit, which is 
minimal at optimal doping but increases up to the normal- 
state value in the strongly under- and overdoped regions of 
the phase diagram. This residual term is likely to have an 
intrinsic origin and is not due to, e.g., residual flux in the 
crystals. The values of 7 r imply that the superconducting 
fraction is a strong function of doping and is only close 
to unity near optimal doping, which is similar to the be- 
havior of cuprates and warrants further theoretical study. 
A detailed analysis of the specific heat shows that the 
superconducting state is described well using a two-band 
model with isotropic s-wave superconducting gaps, which 
are interpreted in terms of a more microscopic three-band 
picture. At optimal doping there is a tendency to strong 
coupling for the larger of the two gaps, whereas the smaller 
gap is nearly independent of doping. 



This work is part of the DFG priority program (SPP 
1458) " High- Temperature Superconductivity in Iron Pnic- 
tides" funded by the Deutsche Forschungsgemeinschaft 
(DFG). 
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Fig. 2: (a) Phase diagram of Ba(Fei_ 2 ,Co :c )2As2 derived from 
our specific-heat (circles) and thermal-expansion (squares) 
measurements, (b) Sommerfeld coefficient 7„ and the resid- 
ual density of states 7r. Measurements of Mu et al. [28] are 
shown for comparison, (c) Normalized specific-heat jump. The 
dashed line represents the BCS value AC/7 n T c =1.43. (d)-(e) 
Reduced gaps and band densities of states (normalized to 1 
mole of superconducting material) derived from the two-band 
analysis. Solid lines are guides to the eye. 
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Doping evolution of the superconducting gaps and electronic densities of states in Ba(Fci_ :E CO; E )2As2 
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Fig. 3: Evolution of the superconducting gaps with critical tem- 
perature T c for Ba(Fei_ :E Co :E )2As2. Filled and empty symbols 
represent overdoped and underdoped samples, respectively. 
Gap amplitudes of Bao.6sKo.32Fe2As2 [Tl] are also shown (blue 
symbols). Lines are guide to the eye. 
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